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1. Introduction 


The statistician concerned with economic 
entomology is frequently asked to recommend 
a sampling technique for estimating the den- 
sity of an insect population. The principles to 
be adopted in devising a suitable technique 
are now fairly well understood; their object is 
to obtain, by a method economical of time and 
effort, an estimate known to be sufficiently ac- 
curate for practical use, and, as a secondary 
consideration, some measure of that accuracy. 
(The intervention of the statistician has some- 
times resulted in greater attention being given 
to standard errors than to the estimate of popu- 
lation!) The difficulties are usually more en- 
tomological than statistical, for the technique 
can only be decided in relation to the biology 
and ecology of the insect studied. Different 
species of the same genus, or even different 
stages in the life cycle of a single species may 
require very different sampling methods. 

In this paper, attention will be restricted to 
a consideration of statistical problems arising 
in the development of a sampling technique for 
soil insects, and, as an illustration of these, the 
estimation of wireworm populations will be dis- 
cussed in some detail. 


2. Principles of Sampling 
Even within a single field which provides 


conditions suitable for the life of a particular 
species of soil insect, the population density 
may vary widely from point to point. At any 
one point, the density may change with time, 
showing cyclic and secular trends as well as 
sudden changes induced by ploughing or 
other violent alterations of environment. Con- 
sequently, unless the sampling procedure in- 
volves counts of randomly or, at least, objec- 
tively selected portions of the population, the 
estimate may be seriously biased. 

The sample used for estimating the mean 
population density in a field generally com- 
prises a set of sampling units, small circles or 
squares of surface dug to a depth great enough 
to include all, or nearly all, insects of the 
species studied which are within the boundary 
of the unit. The estimate of population is 
then calculated from the mean number of 
insects per unit. The method of extracting 
insects from the soil, whether it be hand-sorting 
or a mechanical process, should be capable of 
discovering almost all of the species which 
are present in the sample, as inefficiency in 
extraction may vary from field to field, from 
day to day, or from worker to worker, thus 
ruining the comparability of results. Herein 
lies the chief criticism of any attempt to esti- 
mate populations from baiting records: the 
area within which the attraction of a bait will 


operate cannot be exactly defined, and there 
can be no certainty that all insects in the 
vicinity will respond to the attraction. 


When a sampling procedure is being de- 
veloped, its precision must be examined; 
sampling units must then be taken at random, 
either from the whole field or in equal num- 
bers from every one of a number of equal areas 
into which the field has been divided. The 
standard error of sampling is estimated from 
discrepancies between counts in units from the 
same part of the field. After the establishment 
of a sufficiently accurate technique, there is no 
need to calculate errors for every field sampled, 
but an occasional check should be made; in- 
deed, the precision will be increased, though 
it can no longer be measured, by dividing the 
field into as many areas as sampling units and 
taking one unit from a random point in each 
area (5). 

Clearly one important consideration is the 
size of the sampling unit. The labour of 
extracting the insects from the soil will be 
roughly proportional to the volume of soil, or, 
since the units are dug to a fixed depth, to the 
total surface area of the sample. Also, if 
samples have to be taken from field to lab- 
oratory for examination, this total area will be 
an important factor in the planning of trans- 
port. To dig many small units probably re- 
quires longer time than to dig a few large ones 
which cover the same total area, but, as a first 
approximation, the efficiencies of units of dif- 
ferent sizes may be compared in terms of the 
total sample areas needed to give equal pre- 
cision in the population estimates. For ex- 
ample, if units with surface areas wi, ws have 
sampling variances per unit of 0, vs respec- 
tively (measured on the same scale of popu- 
lation density), samples of total areas A, As 
will lead to estimates of equal precision if 

U1 vi/ Ax = usv2/ Ao, 
and the efficiency of the second unit relative 
to the first is then Ai/Ao. 

If individual insects are distributed in the 
field entirely at random and independently of 
one another, the number per sampling unit 


will follow a Poisson distribution. The vari- 
ance, o*, of the number per unit (of area uw) is 
then equal to the mean number, 4; an easily 
verified consequence is that the variance of an 
estimate of population based on 7 units is 
inversely proportional to nu, the total area of 
the sample, so that all sizes of unit are equal in 
efficiency. Any departure from complete ran- 
domness will usually be in the direction of 
making o* greater than w— though possibly 
some form of competition might lead to a more 
even distribution than the Poisson, thus redu- 
cing o°— and the smallest units are then al- 
ways the most efficient. The magnitude of the 
differences in efficiency required to be investi- 
gated, however, in order to assess how far the 
increased efficiency of a small unit compensates 
for the greater labour in the field. 


3. Sampling for Wireworms* 


In recent years much attention has been 
given to problems connected with the esti- 
mation of wireworm populations, and the 
development of this work provides a good 
illustration of the use of sampling methods for 
soil insects. The first use of a soil sampling 
method for wireworms was probably that of 
Roebuck (4), but he did not investigate the 
precision of the method. Jones (1) reported 
the results of sampling about 50 fields, using 
sampling units of 1, %4, or 1/16 square foot 
and 25, 50, or 100 units per field; he consid- 
ered that his counts agreed with the Poisson 
law. On a priori grounds this conclusion 
seems unlikely to be correct, for the younger 
larvae might be expected not to have dispersed 
completely from points of oviposition and vari- 
ations in soil characteristics and plant cover 
might also interfere with a completely random 
distribution. Cochran’s examination of samp- 
ling results obtained by Ladell (2) verified 
that the error variance was generally greater 
than the mean count, and, as will be seen 
below, more extensive work has since con- 
firmed this. 

When the campaign for the ploughing of 
old grassland began in 1939, it was realized 
that the succeeding arable crops were likely 


*The wireworm is the larval stage of the elaterid or ‘click’ beetle. The most important British 
species are Agriotes lineatus, A. obscurus, and A. sputator, though other genera are also found, 


notably Athous and Cryptohypnus. The 
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make wireworms a major pest of crops grown in the few 


to suffer severe wireworm damage. The Ad- 
visory Entomologists of England and Wales, in 
the thirteen Provinces into which the country 
is divided for agricultural advisory work, 
therefore began a vast program of sampling, 
the Wireworm Survey of England and Wales 
(3). The purpose of this Survey was to dis- 
criminate between those fields in which careful 
choice of crop and cultivations would be need- 
ed in order to combat the wireworm danger 
and those in which low populations were un- 
likely to harm any but the most susceptible 
crops, and also to obtain further information 
about wireworm populations and their effect on 
crops. At first the standard technique was to 
take twenty 6 in. square units from a field of 
10 to 20 acres, selecting a random pair of 
points within each one-tenth of the field area. 
Later, in order to make possible the sampling 
of a larger number of fields, the amount of 
soil per field had to be reduced. As a cylind- 
rical core of 4 in. diameter ( 1/500,000 acre) 
had by then been shown to be more efficient 
than the 6 in. square, a standard of twenty 4 
in. cores was used; rather more cores were 
taken from larger fields, approximately in 
proportion to the square root of the area, since 
the importance of accurate classification in- 
creased with the size of the field. The units 
were dug to a depth of 6 in., and twenty of the 
standard cores weighed about 1 cwt. In gen- 
eral the larvae were extracted from the soil by 
hand-sorting; only those exceeding 5 mm. in 
length were scored, as many of the smaller 
ones were likely to escape notice. In this way 
25,000 fields were sampled before the end of 
1942, and since then the total has probably 
increased to over 40,000. The figures dis- 
cussed here refer only to larvae of species of 
the genus Agriotes, which formed by far the 
Major part of the population in most of the 
fields. 

After the accumulation of sufficient evidence 
on the sampling errors, the requirement of 
random selection of sampling points was 
dropped. Some Advisers then adopted sys- 
tematic patterns of points, but no attempt at 


standardization was made and patterns were: 


varied to suit the shapes of fields. A common 
scheme was to sample in two lines parallel to 
the sides of the field and as far apart as was 
each from the nearer side; ten sampling 


points were equally spaced along each line, or 
displaced laterally a fixed distance alternately 
to either side of these positions, These sys- 
tematic patterns might lead to biased estimates, 
but the earlier work had shown no evidence 
that, for example, the population density at 
the edges of a field differed consistently from 
that at the centre. A slight bias was unimport- 
ant for the purposes for which the estimates 
were to be used, since all fields would be af- 
fected to a similar degree. No attempt was 
made to compare the merits of different sys- 
tematic patterns, but those which are well 
spread over the whole field are unlikely to dif- 
fer in precision markedly and consistently 
from field to field. 


4. Precision of the Estimates 


For the first fields sampled in the Wireworm 
Survey, the values of s?/m were calculated 
(m is the mean count per sampling unit and 
s* the observed variance between sampling 
units) ; in general this quantity exceeded un- 
ity, and tended to increase with increasing m, 
thus indicating that the distribution of larvae 
within a field did not satisfy the Poisson law. 
An alternative test of the Poisson distribution, 
which was found more convenient, is to plot 
s/m (the coefficient of variation) against m; 


the ordinate should average 1/\/m if the distri- 
bution law is the Poisson, but if the larvae are 
not distributed entirely at random the ordinate 
will be greater. 

In order to estimate the average precision of 
sampling at different levels of population, and 
to assess the relative merits of different samp- 
ling units, the relationship between s and m 
had to be studied. It must be remembered 
that the problem of advising farmers on war- 
time management and cropping was urgent, so 
that neither time nor staff could be spared for 
a preliminary extensive investigation of samp- 
ling variation. The Poisson distribution was 
soon found to give too low a variance (5); an 
empirical relationship had therefore to be 
sought, and an intelligent guess at this had to 
be based on the results of the first few (about 
100) fields sampled. Bartlett (6) has pro- 
posed that for data of this type the sampling 
variance may be related to the mean by the 
equation 

s* = am + bm’, 


and Bliss (7) has proposed the alternative 
3? = am?; 

a later examination has suggested that the 
second of these would satisfactorily fit the 
wireworm data, but, since s/m had been plot- 
ted against m in testing the Poisson hypothesis, 
the relationship was originally graduated by 
drawing a freehand curve on this diagram. 
Fields having the same value of m often showed 
wide variations in s, so that the exact form of 
the relationship could not be very accurately 
determined nor would it, if known, be a reli- 
able guide to the sampling errors in individual 
fields. Nevertheless, as further data were col- 
lected, the original freehand curve was found 
to need only slight modification. 

In the winter of 1940-41, 2272 fields which 
had been under grass in 1940 were sampled. 
From these the relationship between s and m 
could be estimated much more satisfactorily; 
for practical purposes, only a smoothed em- 
pirical curve was required and the true alge- 
braic form was of no particular interest. 
Values of s/m corresponding to a given m 
were averaged and plotted against m, and 
again a freehand curve was drawn through the 
points; this curve (Fig. 1) shows the sampling 
errors to be substantially greater than for a 
Poisson distribution except at very low values 
of m. Below m=0.8 the ordinate is approxi- 
mately (1-++0.28m)/\m, above m=1.0 it is 
approximately 0.60+0.62/m, and these form- 
ulae were used for interpolation. 

TABLE 1 
Efficiencies of Various Sampling Units 
Relative to the 4 in. Core 
Population density Per cent efficiency of 


No. per 20 1,000 6 in. 2¥2 in. 2 in. 
4 in. cores per acre square core core 
1 ao 92 100 100 

2 50 75 101 101 

4 100 69 102 102 

12 300 65 106 107 
20 500 62 lll 113 
40 1,000 59 120 126 
80 2,000 58 131 139 


Results for other years agreed closely with 
those shown in Fig. 1. Similar curves have 
been obtained for other sizes of sampling unit, 
and from them have been estimated the aver- 
age relative efficiencies of the units at different 


levels of population (3). For example, for a 
population of one million per acre, Fig. 1 shows 
s/m for 4 in. cores to be 0.91, and the corres- 
ponding value for 2 in. cores is 1.62; the larger 
core has four times the area of the smaller, 
and the efficiency of the larger relative to the 
smaller is therefore 4x(0.91/1.62)? = 1.26. At 
low populations, since the Poisson Jaw is then 
almost satisfied, all units are of about equal 
efficiency, but, as Table 1 shows, at high popu- 
lations the differences become important. The 
4 in. core is much to be preferred to the 6 in. 
square, and the 2% in. core has definite advan- 
tages over the 4 in., though two and a half 
times as many cores are required to give the 
same total area. The 2 in. core is very little 
superior to the 2% in. and the difference 
scarcely repays the trouble of taking one and 
a half times as many cores. The ideal is pos- 
sibly somewhere between 4 in. and 2% in., the 
soil type being a factor affecting the choice, 
since small cores are difficult to use on heavy 
or stony soils. It may be noted that Ladell’s 
results and, on re-examination, Jones’s also 
are in agreement with these from the Wire- 
worm Survey. 

In the winter of 1941-42, efforts were made 
to obtain estimates of sampling variation from 
fields still in stubble after a cereal crop in 
1941, as it was suspected that errors of samp- 
ling might be larger under those conditions on 
account of a tendency for the larvae to con- 
gregate in the rows. Results from 262 fields 
as is seen from the curve relating s/m to m, 
also shown in Fig. 1. At low populations, 
sampling in stubble is as accurate as in grass, 
but at populations of about 500,000 per acre 
15 per cent more cores are needed in stubble 
to give the same accuracy as in grass; the 
necessary increase in cores rises to 50 per cent 
at populations of one million per acre, and 
may be as much as 100 per cent at two million. 

When the curves in Fig. 1, and others simi- 
lar, had been established, they were used to 
give the standard error and fiducial limits of 
any estimate of population (3,5), instead of 
assigning to each estimate an error calculated 
from that sample alone. Thus any real dif- 
ferences in sampling variation between fields of 
the same average population were ignored. In 
the main, advice on the cropping of a field had 
to be based on the estimated population, and, 
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STANDARD ERROR PER 4"CORE AS PROPORTION OF MEAN POPULATION DENSITY 


10 20 30 49 


MEAN POPULATION DENSITY 
Fig. 1: Relationship between sampling error 
and mean population density for wire- 
worm sampling by 4 in. core. 
*means for 2272 fields under grass in 1940 
and sampled in 1940-41. 


providing that the technique and intensity of 
sampling were such as had previously been 
shown to give satisfactory accuracy on the 
average, the standard error of the individual 
estimates was not of great importance. When 
errors did not have to be calculated for every 
field, systematic sampling plans could be adopt- 
ed with reasonable assurance that the errors 
would not be greater than predicted from Fig. 
1, and all units from one field could be bulked 
if ease of transport and examination made this 
desirable. Bulking did, indeed, destroy the 
possibility of discovering fields which were 
heavily infested at one end, lightly at the 
other, but previous experience had shown 
that a difference of this kind sufficiently great 
as to justify different advice being given for 
the two parts occurred very rarely. 
5. Use of the Estimates 

The uses made of the estimates of population 

obtained in the Wireworm Survey were of two 


main types. Only the briefest of comments 
can be made on each type here, but a full 
account has been given elsewhere (3). 

Firstly the estimates were used for the fur- 
ther study of wireworm populations per se, 
and of the factors influencing them. Geo- 
graphical trends in population density and 
associations between soil characteristics and 
population were investigated, as also was the 
effect of arable cropping in reducing the popu- 
lation. Table 2 shows that even high popu- 
lations are reduced to comparatively low aver- 
age figures by two successive arable crops, and 
more extensive data from fields sampled in 
two consecutive years only have been used to 
construct curves which give, for any initial pop- 
ulation, the expected population after one year 
of cropping. Comparison of results for indi- 
vidual fields with the expectations derived 
from these curves can then be used to indicate 
whether particular crops or cultivations are 
more or less effective than average in reducing 
the population. 

TABLE 2 
Mean Populations in Three Successive Winters 
of 104 Fields Ploughed from Grass in 1939-40 
and thereafter under Arable Crops 
Mean Population (1,000 per acre) 


No. of fields 1939-40 1940-41 1941.42 
52 160 170 100 
23 440 360 170 
15 810 440 230 
14 1,380 710 270 


Secondly the relationship between popula- 
tion density and the success of crops was 
studied. In general, detailed developmental 
study of a crop and accurate estimation of its 
yield were impossible, but information on 
plant density and yield samples were obtained 
from some fields. Efforts were made, however, 
to have as many crops as possible graded at 
harvest into one of three grades, and from 
these gradings tables such as Table 3 were 
compiled. These tables enabled the resistant 
and the susceptible of the more common crops 
to be discriminated very satisfactorily. 


Application for entry _as second-class matter is pending. The Biometrics Bulletin is pub- 


lished six times a year—in February, April, June, August, October and December—by the Amer- 
ican Statistical Association for its Biometrics Section. Editorial Office: 1603 K Street, N. W., 
Washington 6, D. C. 

Membership dues in the American Statistical Association are $5.00 a year, of which $3.00 
is for a year’s subscription to the Quarterly Journal, fifty cents is for a year’s subscription to 
the ASA Bulletin and members who pay $1.00 additional receive a year’s subscription to the 
Biometrics Bulletin. Dues for Associate members of the Biometrics Section are $2.00 a year, 
of which $1.00 is for a year’s subscription to the Biometrics Bulletin. Single copies of the 
Biometrics Bulletin are 60 cents each and annual subscriptions are $2.00. Subscriptions and 
applications for membership should be sent to the American Statistical Association, 1603 K 
Street, N.W. Washington 6, D. C. 


TABLE 3 
Classification of 1,513 Crops of Spring Oats 
in 1941 on Fields under Grass in 1940 


Wireworm population 
1940-41 (1,000 per acre) 
- 300 


889 698 

325 - 600 401 256 

625 - 1,000 176 90 

1,025 - 47 17 
6. A Warning 


For simplicity of presentation, the above 
account of the Wireworm Survey has been 
written as though the development of the 
Sampling Technique preceded and was en- 
tirely separate from its application. In fact, 
the urgent need of giving farmers practical 
assistance in combatting the wireworm re- 
quired that the study of sampling, the investi- 
gation of relationships between wireworms, 
soils, management, and crops, and farm advis- 
ory work should proceed simultaneously. This 
had always to be borne in mind when attempt- 
ing to use the records to give information of 
research value. The fields sampled were not 
’ randomly selected, but were chosen in response 
to requests for advice, and comparisons be- 
tween mean populations for different regions 
or soil types may therefore be biased, though 
qualitatively if not quantitatively the main 
conclusions reached are likely to be correct. 
Again, records of fields where advice had been 


Number of crops graded 
Total Satisfactory 


Percent 
Poor Failed unsatisfactory 
156 35 21 
115 30 36 
55 OL 49 
17 13 64 


given in one year were for improving the next 
year’s advice. An interesting consequence of 
the type of advice given is that the effect of 
high populations, as shown by tables such as 
Table 3, may apparently diminish from year to 
year. For example, at high populations wheat 
would be considered undesirable for poor 
farmers or on poor land, but good farmers 
would be advised that with careful-cultivation 
and generous seeding they had a good chance 
of a successful crop. Consequently comparison 
of wheat crops on land at high and on land at 
low levels of population would tend to under- 
estimate the true effect of wireworms in re- 
ducing the crop. As the Survey developed, 
such advice would be given more confidently 
and farmers would be increasingly willing to 
accept it, so that the apparent difference in 
results for wheat at extreme levels of popu- 
lation might be expected to diminish. Table 
4 gives some evidence of the occurrence of 
this phenomenon. 


TABLE 4 
Percentages of Unsatisfactory Crops of Winter Wheat 
amongst All Crops Graded in 1941, 1942 and 1943 


Wireworm population 


Percentage unsatisfactory 


(1,000 per acre) 1941 1942 1943 
- 300 28 yaa 21 
BPAY Te 600 42 31 25 
625 - 1,000 50 46 32 
1,025 - 61 48 40 
7. Summary 


The general problem of taking soil samples 
for estimating the population density of a 
soil insect has been considered. As an illus- 
tration, the technique used in the Wireworm 
Survey of England and Wales has been dis- 
cussed. Though the figures analysed refer to 
one genus only, and exclude wireworms less 
than 5 mm. in length, the main findings are 
probably of wider application. Sample counts 
show the distribution of wireworms in a field 
to be not entirely random, so that small samp- 
ling units are more efficient, volume for volume, 
than large. The uses made of the estimates of 


wireworm population have been briefly indi- 
cated, and a caution about the uncritical draw- 
ing of conclusions from surveys of this type 
has been given. 

The data on which this paper has been based 
were obtained by the Advisory Entomologists 
of England and Wales; most of the numerical 
results have already been published in refer- 
ences (3) and (5), but a small amount of 
previously unpublished material has been in- 
corporated. I should like to express my 
thanks to all concerned for making these data 
available to me. 
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THE ESTIMATION OF VARIANCE COMPONENTS 
IN ANALYSIS OF VARIANCE! 
S. LEE Crump 
Iowa State College, Statistical Laboratory 


1. Introduction. Since its introduction in 
1925 by R. A. Fisher (7), the analysis of 
variance has been used most widely to obtain 
tests of the significance of treatment effects. 
When he introduced this technique Fisher 
(7 § 40) indicated a further use of the analysis 
of variance. If an observed statistical vari- 
ate, e.g. the plot yield of a varietal experi- 
ment, is assumed to be the sum of several 
separate effects (variety, block, etc. in the 
case of the varietal experiment), the variance 
of each effect will contribute to the total 
variance. The second use of the analysis of 
variance provides estimates of these several 
variance components. It is the purpose of 
this discussion to point out the hypotheses 
appropriate to the two uses of the analysis 
of variance and to explain its use to estimate 
variance components. 

2. The Hypotheses. Although the hypothe- 
ses in each case are based on the same funda- 
mental equation, they are essentially quite 
different. In order to illustrate, let us con- 
sider a randomized blocks experiment with 
treatments Ai, As, --+, Ag arranged in blocks 
Bi, Bz, «++, By and with 7 observations taken 
in each treatment in each block. If we denote 
by aij the j** observation in the it block for 
the h** treatment, we have for a fundamental 


equation 

(1) Ynij =btbant+Bit oBnrtFniz- 
a aa 
Oey 
ii =1, 2, hyn 


1Journal Paper No. J-1331 of the Iowa Agricultural Experiment Station, Ames, Iowa. 


In equation (1) » denotes the effect common 
to every observation, i.e. the general mean 
of y, ax the effect common to every obserya- 
tion in Ay, 6; the effect common to every 
observation in B;, a8;; the effect common to 
every observation in both A; and B,, and ji; 
the random effect peculiar to the hij** observa- 
tion. 

In the first use of the analysis of variance, 
we have the following assumptions: 


i. The ¢j:; are normally and independently 
distributed with mean zero and vari- 
ance o;’, 

ii. The an, Bi, aBn: and w are parameters 
which remain fixed from sample to 
sample. 


The problem here is to estimate the ap, 6:, 
o8,; and w and to test the null hypothesis 
that any set of these parameters, a, a2, -+-, 
@a, Say, are all equal to zero. The analysis 
of variance provides a solution to this prob- 
lem. 

In the second use of the analysis of vari- 
ance, the assumptions are as follows: 


i. The aa, Bi, aBar and pi; are all ran- 
dom variables independently distributed 
about means zero. 

ii. The parameters in this case are the 
variances of an, Bi, oByi and Shih which 
we denote by a7, og”, o°xg and ox? respec- 
tively, and pz. 


Here the problem is to estimate og2, 0%, ag 


Project No. 890, 
a 


TABLE 1 
The Analysis of Variance for a Randomized Blocks Experiment 


Degrees of Mean Average value of 
Source of variation freedom square the mean square 
Treatments (A Classes) a—1 M.S.4 of +n egt nba 
Blocks (B Classes) b-1 M.S.2 o?+No79g+naag* 
AB Interaction (a—1)(6—1) M.S.48 oP + Nere8 
Error (within Subclasses) ab(n—1) M.S.z 3? 


and o;*. It should be noted that no assump- 
tion is made about the form of the distribu- 
tions. . 

3. The Average Values of the Mean Squares. 
Table 1 shows the analysis of variance for our 
randomized blocks example. Under the as- 
sumptions of the second hypothesis above, it 
may be shown that the average values, in 
repeated samples, of the mean squares in 
column 3 are the expressions shown in the 
fourth column of this table. 

A consideration of these expressions sug- 
gests that we may be able to formulate a 
general rule for determining the average 
values of any mean square from a multiple 
classification with equal subclass numbers. 
We will give such a rule. Consider an 
ABCD - classification with a A-classes, 
b B-classes, etc. and with » observations in 
each of the smallest subclasses. As in our 
preceding example, we write for a funda- 
mental equation 


Vrijkee =ptantBit *+++aBiiteaynjt:*: 
+ oBynigtoaBoniet +++ oP yank. + onijk-- 
B= 2 eee 
i PP RCS 
etc. 


where a, is common to all individuals in the 
h* A class etc. As before let the variances 
of an, Bi, aBri, «** be oo’, of*, cag, :*-- In 
the analysis of variance for this classification, 
denote the mean squares in a manner corre- 
sponding to that of Table 1. With this 
notation the rule may be stated: ‘‘The average 
value of M.S.z (mean square within sub- 
classes) is oy*. The average value of any 
other mean square is a linear combination of 
oy? and all other o?’s whose subscripts contain 
every greek letter corresponding to the sub- 
script of the mean square in question. The 
coefficient of o;? is unity. The coefficient of 
any other o? is A=nabcd---- divided by the 
small roman letters corresponding to the 


greek subscript of the o%.’”’ Let us illustrate 
the rule with two mean squares from a four- 
fold ABCD classification in which a=2, b=3, 
c=3,d=4andn=2. WehaveA=2-3-3-4-2 
=144. According to the rule the average 
value of M.S.4z¢ contains o;?, ogy and og8y5- 
Thus, we have 


E(M.S.azc) =e nr ha Fi 


144 
meds nee 3.37 173. g.47 2878 


= 07+ 80%x8y~t20%a875 

where E( ) denotes average value of the 
term in the parentheses. Similarly 
E(M S.4c) =072+2402xy+80%a8y 

= 60728 ae 207 aBy6- 
This rule will apply to any multiple classi- 
fication with equal subclass numbers. In the 
case that =1 the rule is applied with a; set 
equal to zero. In this case o°ggs... assumes 
the role of o;?. 

4. Estimation of the Variance Components. 
Returning to column 4 of Table 1 we see that 
E(M.S.4s—M.S.z) =n07a8, 
E(M.S.3—M.S.48) =naeg’, 
E(M.S.4—M.S.48) =nbo,?2, 
and hence that unbiased estimates of the o*’s 
may be obtained from linear equations in the 
mean squares. If we let the “hat” (“4) de- 
note “estimate of,’’ we have for example 


(2) 


é4eg=2(M.S.ss—MS.2), 


Gt = (M.S.4—M Sis). 


In any multiple classification, when the aver- 
age values of the mean squares have been 
written down, the linear combination of the 
mean squares which estimates a given o* or 
variance component will be evident. 

In most cases it is desirable to know the 


variance of an estimated component. If each 
of the random elements in our fundamental 
equation follows a normal distribution, it may 
be shown that any analysis of variance mean 
uh 
average value of the mean square in question, 
f the corresponding degrees of freedom, and 
x° follows the ordinary type III distribution 
with f degrees of freedom. Hence, since the 
variance of x? is 2f, the variance of any mean 


square is distributed as where oo? is the 


. 2oo4 a 
square is =e Further, it may be shown 


that the mean squares are independently dis- 
tributed, and thus we may write the variance 
of any linear function of the mean squares. 
For example consider the first of equations (2). 
Let 

E(M.S.4z2) = of + Noeg = 00", 
and 


— or. 


E(M.S.z) =o? 
Then since 


aS *(MS.sp— M.S.z) 
we may write 


A 2 oot ot 
Vers) =5(F+), 

where fo=(a—1)(6—1) and fi=ab(n—1), 
Now o? and o are unknown population 
values and we will obtain biased estimates of 
V(x) if we substitute M.S.4z and M.S.z 
for oo? and oa}? respectively. To correct for 
this bias, Daniels (4) gives the following for- 
mula to estimate V. 


Sp seed (MSPs S22) | 
GHGs) =3( at te =) 


The variance of estimates of other compo- 
nents may be estimated in a similar manner. 
In large samples the estimates are normally 


+ 


distributed under the condition of the pre- 
ceding paragraph. Thus fiducial limits may 
be placed on the estimates when large samples 
are available. 

5. An Example. We shall now illustrate 
the procedures described in the preceding 
sections with a numerical example. Our data 
are drawn from a series of genetic experiments 
on egg production and comprise the total 
number of eggs laid by each of 12 females 
from 25 races of Drosophila melanogaster, on 
the fourth day of laying, the whole experi- 
ment being carried out 4 times.2 The anal- 
ysis of variance for this set of data is shown 
in Table 2. 

The rule for writing the average values of 
the mean squares given in section 3 may be 
verified for the last column of Table 2. Fol- 
lowing equations (2) we obtain the following 
estimates of the variance components: 


$2 =M.S.z=231, 


» 


A 1 
oY¢p = 79 (M-S.ze>M.S.z) = 19, 


3,2 = AMS.2—-MS.rp) = 58send 
Gi eSiy = Mi SIa) S152 


~ 300 
Applying the principles used in obtaining 
equation (3), we have the following estimates 
of the variances of our estimated components: 


Voz) = 22) i sa =e) 
ia ia lear HD) =e; 
Pe) maa (Mag ta) = 354, 
VG) =55 3 ae + Mian) ae 


?T am indebted to Dr. J. W. Gowen of Iowa State College for permission to use these data, 


TABLE 2 


Analysis of Variance of Total Egg Production of 12 Females 
(D. melanogaster) from 25 Races in 4 Experiments 


Source of Degrees of Average value of 
variation freedom Mean square the mean square 
Experiments 3 M.S.2 =46,659 a+ 1207,,+30002 
Races 24 M.S:n=0 3,245 a+ 120?,, +40," 

EXR 72 M.S.cr= 459 o¢’+ 1207. 
Within Subclasses 1100 M.S.z = 231 oe 


One may ask what sort of question may be 
answered with the help of the estimated 
variance components. We shall illustrate 
with one question. Suppose that it is desired 
to estimate the mean egg production of the 
a race on the fourth day of laying to some 

_ specified degree of accuracy. Let %.;. be the 
mean egg production estimated from n fe- 


males in each of e experiments. Then with 
the notation of section 2, we have 
Pr mph ee Oe oh ten 


Sia: : ‘rhein 
en 


- 


Hence we have for the variance of #.;. about 
its average value, utp; 


1 1 
V@.:.) = (oe +o%ep) +7 (er) 
and we may estimate this quantity by 


“ 1 1 
(4) V(E.i-) a 19) +5, (230). 
From equation (4) it is clear that, by increas- 
ing the number of experiments sufficiently, 
we may reduce V(Z.;.) to any desired level. 
However, increasing the number of females 


per experiment indefinitely still leaves us 
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with V(z.;.) =e In practice, the decision 


to increase 7 or e or both enough to reduce 
V(é.:.) to the desired level will depend on 
the relative costs of these alternatives. 

6. Estimation of Variance Components in 
More Complicated Cases. When the classes 
have unequal numbers of individuals, it be- 
comes extremely tedious to work out the 
average values of the mean squares. The 
principles, however, remain the same, involv- 
ing straight forward algebra and elementary 
probability laws. Winsor and Clarke (14) 
have given the results for a one-fold classi- 
fication into groups with unequal numbers. 
Consider data arranged in a groups with n; 
observations in the i** group (¢=1, 2, «+-, a). 


a 
Let N= D7n;. In our previous notation 
i-1 
M.S.4 and M.S.z denote the mean squares 
between and within groups respectively. 
Winsor and Clarke give the following average 


values for M.S.4 and M.S.z: 


E(M.S.z)=o7 and E(M.S.4)=¢;?+m000%, 
where 


1 ws 
mom (WN —y2") 

The average values of the mean squares for 
an extended ‘‘groups within groups” type of 
classification have been given by Ganguli (8), 
Hetzer, Dickerson and Zeller (9) and Finkner, 
Morgan and Monroe (6). 

When the classification is of the ABCD ---- 
type with unequal subclass numbers, there 
are a number of methods of analysis available, 
e.g. the method of fitting constants, Yates 
(15), the method of weighted squares of 
means, Yates (15), and the method of ex- 
pected subclass numbers, Snedecor (12). 
Each of these methods of analysis can give 
unbiased estimates of the variance com- 
ponents. At the present time, however, it 
is not known which method gives ‘‘best”’ 
estimates. 

7. Other Literature on the Estimation of Var- 
tance Components. There are certain rela- 
tionships between the results of the analysis 
of variance under the two hypotheses of sec- 
tion 2 which sometimes enable us to interpret 
tests of significance in terms of variance com- 
ponents. Discussions of such interpretations 
are given by Cochran (1), Crowther and 
Cochran (3), Wilm (13) and Yates and 
Cochran (17). 


The use of the analysis of variance to esti- 
mate variance components has wide applica- 
tion in the selection of efficient sampling 
designs. Excellent discussions of these ap- 
plications are found in Cochran (2), Yates 
and Zacopanay (16) and Youden and Mehlich 
(18). 

Another field in which estimates of variance 
components are widely used is Genetics. 
Dickerson (5), Hetzer, Dickerson and Zeller 
(9), Lush and Molln (10) and Sprague (11) 
provide illustrations of applications in this 
field. 

8. A Note of Warning. It must be re- 
membered that in using the analysis of vari- 


ance to estimate variance components, we — 


have assumed the elements of the funda- 
mental equation to be randomly selected from 
an infinite population. In an experiment 


10 


a 


where three widths of spacing some crop are 
purposely selected for trial, it is not reason- 
able to regard these widths as random samples 
from all possible widths. On the other hand 
the blocks in a field experiment may some- 
times quite reasonably be regarded as a ran- 
dom sample of all such blocks. In sampling 
production from say three machines in a 
factory, where these machines constitute all 
the machines which the factory has or is 
likely to have, it is more reasonable to regard 
these machines as the whole of a finite popu- 


lation than to consider them as random 
samples from some infinite population. If 
the factory owner is sampling production with 
a view to purchasing more machines of the 
same type, the three machines may be appro- 
priately regarded as samples of the infinite 
population made up of all machines of the 
same type. Daniels (4) treats the case of 
finite populations. The assumptions must 
be considered carefully in any problem where 
it is desired to estimate variance components. 
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ANALYSIS OF SCORES FROM SMELLING TESTS* 
Witiram DowE.L. BaTEN 
Professor of Mathematics, Michigan State College, East Lansing, Michigan 


Three small bottles containing 1%, 10% and 
30% of rose extract, colored so that the 
strengths could not be detected by sight, were 
given to individuals, who were requested to 
arrange them in ascending order according to 
strength of aroma. To be able to determine 
whether or not those taking the test really 
showed ability to detect the different strengths 
of aroma, it is necessary to determine the num- 
ber of correct arrangements which may be 
made by mere chance. A person, with no 
ability to smell, might by chance alone, arrange 

the bottles in the correct order. 


known how many assortments on the average 
can be made by chance, then it will be possible 
to ascertain whether or not the judges were 
able to distinguish the various strengths of 
this material. 

If the following is the correct arrangements 
of the bottles: 


A B Cc 
1% 10% 30% 


then the following six arrangements, desig- 


nated as (a), (b), (c), (d), (e) and (f), are 


After it is possible: 


Journal article No. 785 (N.S.) of the Michigan Agricultural Experiment Station. 


Ut 


TABLE I 
1% 10% 30% Value 
(a) IG B A 0 
aerate A B 
(c) B Cc A . 
(4) A Cc B 
(cyeaB A Cc 2 
(ye W B Cc 3 


The order (a) means that the judge thought 
that the 1% solution was stronger than the 
30% solution. This is the poorest judgment 
that he could make, although the second 
bottle is placed correctly. Order (b) means 
that the individual thought the strongest was 
the weakest, that the weakest was the next to 
the weakest, and that the next to the weakest 
was the strongest. This decision does not 
appear as bad a judgment as the above where 
the strongest and weakest solutions were inter- 
changed. Order (c) indicates that the samp- 
ler thought that the weakest was the strongest 
concentration, that the next to the weakest was 
the weakest and that the 30% concentration 
was stronger than the 10%. Arrangement (d) 
shows that the smeller became confused on the 
two strongest aromas, but was able to distin- 
guish the weakest from the two strongest. This 
person considered the 10% solution to be 
stronger than the 30% solution. Assortment 
(e) means that there was confusion about the 
correct ranks of the weakest and the next to 
the weakest. In the previous two orders, the 
tester became confused on the two adjacent 
strengths. 

Let a value of 0 be given to (a), a value of 
1 be given to a combination of (b) and (c),a 


x Probability 
0 po” 

1 po” pr 

2 

S 

3n p"s 


value of 2 be given to a combination of (d) 
and (e) and a value of 3 be assigned to order 
(f), with respective probabilities 1/6, 2/6, 
2/6 and 1/6. These weights are based upon 
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the number of interchanges in an assortment 
necessary to secure the correct assortment. 


Let the general case be considered, where 
po, pi, po, and 3 represent respectively the 
probabilities of four events 0, 1, 2, and 3 
occurring (in our case arrangement (a), either 
(b) or (c), either (d) or (e) and arrangement 
(f)). The terms in the expansion of the 
multinomial 


(po+ pit po+ ps)” 
5 


(where the sum takes in all possible occur- 
rences and where mo, m1, “2, and m3 represent 
respectively the number of times each event 
occurred) gives the respective probabilities of 
the various events; that is the first term, fo", 
in this expansion is the probability of event 0 
occurring ” times; the second term, po"1, 
is the probability that event 0 will occur 
nm—1 times and event 1 will occur one time; 


! 
the term, Tisighbioey, is the probability 


that event 0 will occur 7 times, event 1 will 
occur 7 times, event 2 will occur s times and 
event 3 will occur g times. Let the occurrence 
of event 0 correspond to a value of 0, the 
occurrence of event 1 correspond to a value 1, 
the occurrence of event 2 correspond to a 
value of 2 and the occurrence of event 3 
correspond to a value of 3. 
3 
If «= > in;, then x will range from 0 to 3m. 
i=0 

The value x=6, is obtained when m=n—6, 
m=6, m2=73=0, or when m=n—5, m=4, 
m.=1, n3=0, or when m=n—4, m=2, m=2, 
n3=0, or when m=n—3, m=0, m=3, m3=0, 
or when m=n—4, m=3, m=0, m3=1, etc. 
Other values between 0 and 3m can be simi- 
larly analyzed. The following frequency dis- 
tribution gives the various values of x corre- 
sponding to the above events, with their 
respective probabilities or frequencies. 


jbo™bx™pa*pa™ 


no aa i 'n3! 


! 
Gan” eit nb tp: 
! ! 
(n =a 8 jt Gy be tbat mba tbs 


It is necessary to determine the mean 
and standard deviation of this distribution. 
These can be found in the usual way or may 
be secured by taking the first and second 


derivatives, at z=0, of the following expres- 
sion 
Y= (poe™+ pie?+ pre?*+ pre)” 


The first derivative of Y at z=0 is found as 
follows: 


VY'=n (poe*-+ pie?-+ pre? + pre” tat 
X (pie? + 2pre*+ 3pse*), 


V',20=1(Pit2p2+3p3) = p's. 
The second derivative of y at z=0 is 


bo:2=n(n—1)(pit2p2+3p3)? 
+n(pit4p2+9p3). 


The second moment about the mean is 


o22=n[(pit2p2+3p3)(po— p2—2ps) 
+2p2+6p3 ]. 


The standard deviation of the above distribu- 
tion is 


(2) on= Me ee 23) 


and 


+2p.+6p3 ] 


Two hundred people took the smelling test 
described above with the results given in 
Table 1. 

TABLE 1. 


Number of people arranging the bottles with 
arrangements (a), (b) and (c), (d) and (e) 
and (f), where 7 represents the corresponding 
values. 


a Frequency Percentage 
0 18 .090 
1 39 195 
2 73 365 
3 70 .350 


The total of the scores or the value of » is 
395; this total was obtained as follows: 


(OX 18)+(1X39)+(2X73)+ (3X 70) =395. 


Is this total of 395 significantly different from 
the total score which 200 people would have 
made, on the average, if they had no ability 
to smell and had arranged the bottles by mere 
guessing? If chance alone were used, the 
average total according to (1) is 


u',=n(pit2p.+3ps) = 200(2/6+4/6 


since po=1/6, pi1=2/6, p2=2/6 and p:=1/6. 
According to (2), the standard deviation is 


The #-value will show whether or not the 
observed total of 395 is significantly different 
from the average total of 300, which was due 


+3/6) =300, 


to chance. The ¢-value with an infinite num- 
ber of degrees of freedom is 
395 — 300 
pra 035) ae i 


which indicates that the judges were not 
arranging the bottles by chance. ‘The testers 
were evidently using their sense of smell in 
arranging the bottles in the required order. 


Attempts at removing the effects of the last 
aroma. 

Some of the judges said that after they had 
smelled the first bottle their olfactory nerves 
were saturated with the aroma, so that they 
were not able to detect stronger or weaker 
strengths afterwards. To overcome this dif- 
ficulty several materials were used between the 
bottles of rose extract. 

The same test was given to 73 mature 
people with the condition that a bottle of 
peppermint was to be held to the nose before 
smelling the second bottle of rose extract and 
also before smelling the third bottle of rose 
extract. The object of this test was to deter- 
mine whether or not the aroma from the pep- 
permint removed the aroma of the rose extract 
from the olfactory nerves and enabled the 
judge to better distinguish the strength of the 
next bottle of rose extract. In some taste ex- 
periments the tasters eat a cracker, take a 
sip of water, or eat a slice of an apple to re- 
move the effects of the last sample from the 
taste buds. Table 1 contains the expected per- 
centages of the individuals that on the average 
arrange the bottles as (a), (b) and (c), (d) 
and (e), and (f). These percentages were 
obtained from the 200 people taking the test 
and constituted all the information concern- 
ing the ability of people to assort the bottles 
according to strengths of aroma. If the above 
percentages are the true ones, and they are 
used, then the mean and standard deviation of 
the values from 0 to 219, whose frequencies are 
the respective terms of the expansion of the 
multinomial 

(.090 + .195 -+ .365 + .350)™ 
are equal respectively, to 
wx=144 and ox=8.2 

The second column of Table 2 contains the 
results of the test when peppermint was used 
between the bottles of rose extract. The total 
score or value found from these values was as 
follows: 
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(0X9) +(1X13) + (2 X30) + (3X21) =136 


How does this total of 136 compare with the 
expected total of 144 when no peppermint was 
employed? 
The t-value is 
t=—136—144 
8.2 
which indicates that the peppermint did not 


= 098 


better arrange the three strengths of rose 
extract. Some of the judges thought that 
these materials did aid them, others thought 
that they hindered them, and some could not 
detect any assistance or any hindrance. The 
expected frequencies in Table 2 indicate how 
near they are to the observed frequencies; 
clearly the aromas used between the vials did 
not help or hinder the judges in assorting the 
three strengths of aroma. 


TABLE I 


Table 2. Frequently distributions of the values assigned to arrangements (a), (b) and (c), 
(d) and (e), and (f), where i represents the corresponding values 


Materials used between bottles of rose extrace 


Peppermint Turpentine Vanilla Vinegar 
. Observed Exp.|Observed Exp.|Observed Exp.|Observed Exp. 
J freq. freq.|freq. freq.| freq. freq.| freq. freq. 
0 9 6.6 6 6.1 5 6.8 12 5.8 
1 13 142] 11 13 oueds 146 | 10 12.5 
2 30 26.6 | 28 24.8 | 26 27.4 | 23 23.3 
3 21 25.6 | 23 23:8 lana 26.3 | 19 22.3 
No. of people N 73 73.0 | 68 68.0 | 75 75.0 | 64 64.0 
Observed total 136 136 158 113 
Expected total 144 134 148 126 


affect the decisions pertaining to assorting the 
rose extract solutions. 

A bottle of turpentine was used between the 
bottles of rose extract as was done with pepper- 
mint to determine whether or not a whiff of 
turpentine before going to the next bottle in 
the test removed the aroma in the nose. Vanilla 
extract and vinegar were used similarly. The 
results from these tests also are contained in 
Table 2 together with the total scores and the 
expected scores. 

On examining the respective t-values, it is 
found that there were no significant differences 
between the actual totals and the expected 
totals, which show that the turpentine, vanilla 
extract and vinegar did not help the judges to 


It is interesting to note, in the experiment in 
which 200 people were used, that the total 
score or value of 395 was far from a perfect 
score of 600. Thirty-five per cent were able to 
assort all of the bottles correctly; 28.5% as- 
sorted them as either (a), (b) or (c). This 
means that these strengths of 1 per cent, 10 per 
cent and 30 per cent could not be detected ac- 
curately by the majority of the people and that 
people differ considerably in their ability to 
distinguish aromas. 

If a food laboratory desires to select a panel 
of judges for detecting aromas of certain foods 
it would be advisable to carry out a similar 
test and select those who were able to arrange 
the material correctly. 


QUERIES 


QUERY: I have a problem in mind which I 
would like to have answered if possible. 
Source of Degrees of Mean 
variation freedom Square 
Total 99 
Blocks 9 66.667 
Treatments 9 Sept) 
Error 81 1.234 


In the above table of analysis of variance 
both “Blocks” and ‘“‘Treatments”’ are signifi- 
cant at the 1% level. However, the value 
of F for Blocks is twice as great as the value 
of F for Treatments. Does this mean that 
the variance accounted for by Blocks is twice 
that accounted for by Treatments, or might 
the variance accounted for by Blocks and 
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that accounted for by Treatments be approxi- 
mately the same for both Blocks and Treat- 
ments since both are significant at the 1% 
level? 

I am concerned about the interpretation of 
the relative importance of the two significant 
variables. 


ANSWER: Following the ideas developed 
in the paper by this writer in the present issue 
of the Biometrics Bulletin and taking the 
assumptions outlined there as fulfilled, we 
find for estimates of og, o,?, and o2 the fol- 
lowing: 


Ge? = 6.543 
6,2=3.210 and 
¢2=1.234 


where the subscripts 8, 7, and e, refer to the 
blocks, treatments and error components of 
‘variance respectively. Hence, in your ex- 
ample with its small error mean square, the 
estimate of the variance ascribable to blocks 
is slightly more than twice that of the vari- 
ance ascribable to treatments. However, 
these estimates are subject to a large sampling 
error. By approximate methods we may 
place fiducial limits on og? and ¢,2. The 5% 
fiducial limits in this example are 3.423-17.919 
for ég? and 1.650-8.900 for ¢,2.. Any inter- 
pretation of the relative importance of the 
two main effects should take into account 
this large sampling variation. 
S. LEE Crump 


QUERY It is usual to plot the independent 
variable on the X axis and the dependent vari- 
able along the Y axis, but Ezekiel, 1930, 
pp. 129 and 131, reverses the procedure. Since 
T was trying out a linear and a curvilinear 
regression, that change bothered me. Is there 
any rule about this? 


ANSWER It is conventional, as you say, to 
plot the independent variable along the hori- 
zontal axis. Ezekiel apparently does this, the 
amount of water applied being considered the 
independent or controlled variate with cotton 
yield as dependent. 

Perhaps you are thinking that at any partic- 
ular stage of growth the size attained by the 
plant determines the amount of its water re- 
quirement. If so, then you would doubtless 
plot the size of the plants as X along the 
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horizontal axis. Ezekiel’s irrigation problem 
required a different attitude toward the cause 
and effect relation of the variates. 

George W. Snedecor 


QUERY In a covariance problem, the ques- 
tion arose as to how to adjust the treatment 
means for differences in their values of the 
independent variate. Will you explain this 
for a randomized blocks experiment? 


ANSWER Ordinarily the treatment means 
are adjusted by use of the error regression 
coefficient. This results in the set of adjusted 
means, differences among which are tested by 
the usual procedure. The following data are 


needed: 
j. The error regression coefficient, 
b=Sxy/Sx’. 
2. The deviation of the treatment, 


mean of the independent variate 
from the experiment mean. Denote 
this by x. 

. The treatment mean of the depend- 
ent variate designated as Y. The 
adjusted value is then Y—bx. 

In Fisher’s example 46.1, the error row in 
the covariance table has Sx? = 567.5 and 
Sxy = 654.25; hence, b= 654.25/567.5 = 
1.1529. 

If the mean of all values of X is 100 while 
the mean X for a certain treatment is 92.25, 
then x = 92.25 — 100 = — 7.75. 

The mean Y for the same treatment is 82.25. 
Finally, the adjusted mean Y is 82.25 — 
(1.1529) (— 7.75) = 91.18. 

George W. Snedecor 


QUERY In the first query of the August, 
1945, Biometrics Bulletin, page 55, the 
questioner ends his problem with the following 
question: “What . is the probability that 
the . . . samples are not merely samples from 
the same population?” Obviously this ques- 
tion cannot be answered without a@ priori 
knowledge of the probabilities of the popu- 
lations. 

At first thought one might regard the error 
in the question as merely a slip of the tongue 
or poor phraseology. However, recent experi- 
ence with statisticians and experimentalists 
trained in various colleges has convinced me 


that it reflects more often than not a serious 
error in thought. They persistently are found 
writing or saying that the small probability 
derived from a significant test is the proba- 
bility of there being a true difference in the 
populations sampled! 

In fact one need not rely upon inexperienced 
statisticians as a source of such errors. If 
memory serves me right, the same error is 
noticeable in some textbooks. 

The answer to the query gave the probability 
that two independent random samples from the 
same population could produce two inter- 
actions differing (as measured by their ratio) 
more than the ratio obtained from the two 
observed samples. Should you not have re- 
phrased the question so as to make it clear 
what question you were answering? 

ANSWER Yes; and thank you for calling 
attention to the discrepancy. 

Other readers have raised the question as to 
whether I really got at the root of the original 
querist’s difficulty. Of the seven degrees of 
freedom in the 2X2X2 table, only one was 
discussed, leaving three main effects and three 
first order interactions unnoticed. So far as I 
have been able to see, these other six effects 
are not relevant to the query. 

George W. Snedecor 


QUERY Ina cultivation experiment on sugar 
beet, three treatments were tried in seven 
randomized blocks. The analysis of variance 
and the treatment means are as follows: 


Source of Degrees of Mean 
variation freedom square 
Treatment 2 1.534 
Block 6 3.409 
Error 12 0.232 
Mean yields in tons per acre 
Treatment 1 15.88 
2 15.00 
3 15.18 


For treatments, F = 6.61, which does not 
reach the significant value, 6.93. On the con- 
trary, the ratio of the range of means to sx, 
0.88 \/0.232/7=4.84, is beyond its 1% value, 
4.10. This was calculated as indicated in the 
April Bulletin. 

It is my understanding that no significant 
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differences exist if the calculated F value is 
less than that found in the table. Thus, from 
the F test, I conclude that there are no dif- 
ferences, while from the test of range the 
opposite is true. I have heard several explana- 
tions of this discrepancy but would appreciate 
having your comments. 


ANSWER From my point of view, there is 
no discrepancy, because the probabilities 
turned up by the two tests are practically the 
same. If one probability were 0.1 and the 
other 0.01, I should rely upon the probability 
of F as being the more dependable; but I 
should also go over the whole experiment again 
to learn why the two indications were so dif- 
ferent. 

Perhaps your difficulty lies in a somewhat 
arbitrary interpretation of the test. You 
imply that if P is a bit larger than 0.01 there 
are no differences, while if. P is just a little 
smaller than that value, then suddenly there 
are differences. This is not a realistic attitude 
toward a test of significance. There is no 
sharp break in the scale of probability but 
rather a continuous flow. The conventional 
0.05 and 0.01 points are like mileposts on a 
road: you do not jump from 6 miles to 5 as 
you pass the post. Let us review the steps 
taken in a test of significance. 

First, one sets up a null hypothesis in con- 
formity with the objective of his experiment— 
in your case, the hypothesis is that the three 
treatment means are drawn from a common 
population; that is, that the treatments do not 
affect yield. Next, one calculates F for the 
experimental yields. Third, one observes in a 
table the probability of a larger F in random 
sampling from the hypothetical population. In 
your experiment, the probability of a larger F 
is about 0.01; that is, one would expect to get 
a larger F about once per hundred trials if 
the treatments were wholly ineffective. Finally, 
the experimenter makes a decision. He may 
decide that, despite the improbability of his 
sample, he will not reject the hypothesis set 
up: he may not be convinced that treatment 1 
is superior in the population from which the 
experimental sample is drawn. Usually, how- 
ever, with P=0.01, he would reject the null 
hypothesis, concluding that the treatments do 
affect yield, and that in the sampled popula- 


tion treatment 1 is superior to the other two. 

It should be observed that the experimenter 
might make either of these decisions if the 
experiment had resulted in other levels of 
probability. It is customary to decide in ad- 
vance the level that will lead to rejection, but 
it should be clear that the statistical evidence 
is only part (and sometimes a small part) of 
the entire information upon which the decision 
of the experimenter must be based. There is 
no value of P at which it can be said with cer- 
tainty that differences do or do not exist in the 
population. 

George W. Snedecor 


QUERY In a study of the ability of exam- 
iners to predict training success, a two by two 
table was secured for each examiner. Tables 
for eight examiners are given below. The 
problem is to test the hypothesis that the ex- 
aminers do not differ among themselves with 
respect to their ability to predict training 
success, 


Examiner Actual Examiner’s Prediction 
Outcome Pass Fail 

1 Pass 6 2 
Fail 16 20 
2) Pass 11 2 
Fail 19 5 
4) Pass 7 2 
Fail 10 14 
4 Pass fl 1 
Fail 8 a 
59 Pass 3 1 
Fail if ll 
6 Pass 8 0 
Fail 4 12 
7 Pass 5 1 
Fail 9 5 
8 Pass 5 0 
Fail 9 5 

ANSWER Ordinarily to test a 2X2Xr table 


for homogeneity with respect to the r-fold 
classification, one would employ a x? test for 
which 2r—2 degrees of freedom are available. 
In this particular example, however, it seems 
worthwhile to make two tests each involving 
r—1 degrees of freedom as will be explained 
below. 

The sums for the above 2X28 table over 
the 8-fold classification are: 
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Actual Examiners’ Prediction Totals 

Outcome Pass Fail 

Pass 52 9 61 

Fail 82 74 156 
Thus the examiners correctly 


predicted that 52 of the 61 passing students 
would pass, and that 74 of the 156 failing 
students would fail. They obviously did far 
better with the passing students than with the 
failing students. The reason for this is easy 
to surmise—the examiners wished to avoid 
predicting that a passing student would fail 
and in order to avoid doing so, they were will- 
ing to include a number of doubtful cases in 
their passing category. Thus the examiners 
predicted failure only when they felt fairly 
certain that the student would fail; all other 
students were predicted to pass. 

The 9 passing students who were predicted 
to fail were, therefore, gross errors on the part 
of the examiners; while the 82 failing students 
who were predicted to pass merely represent a 
margin of safety used by the examiners in their 
endeavor to include all passing students in 
their passing prediction. 

Two comparisons may be made between the 
examiners. A comparison using only the pass- 
ing students will indicate whether the exam- 
iners differ in their proportions of gross errors 
having selected arbitrarily their own margins 
of safety against gross errors. The x’ test of 
homogeneity on the 28 table involving the 
61 passing students gives x* = 4.7 with 7 de- 
grees of freedom, which corresponds to a prob- 
ability of about .80. 

A comparison using only the failing students 
will indicate whether the examiners differ in 
their margins of safety. Here x°=19.7 with 
7 degrees of freedom and the probability is 
less than .01. Thus while all examiners miss 
a small proportion (about 15%) of the passing 
students, some examiners must pass signifi- 
cantly more failing students than others to 
achieve this small proportion. 

If we may assume that it is easy to predict 
correctly the outcome for very capable stu- 
dents and very dull students, then the exam- 
iners should be judged on their treatment of 
borderline cases. The comparison using onlv 
the failing students enables one to judge the 
examiners on their treatment of borderline 
cases, but this is a circumstance of the data, 


and is possible only because the first x” test 
did not reject the null hypothesis (that there 
was no difference between examiners with re- 
gard to their commission of gross errors). 
If, as might well have happened, some exam- 
iners had used a small margin of safety and 
made few errors in judging the failing students 
while making several errors in judging passing 
students; then it would have been appropriate 
to compare the examiners by the ordinary x* 
test of homogeneity for 2X2r tables. In the 
example at hand, this x? would simply be the 
sum of the two xs obtained above, x°=24.4 
with 14 degrees of freedom, which corres- 
ponds to a probability of about .04. 
A. M. Mood 


QUERY: We are planning to test the adapta- 
bility of some 28 grasses and legumes in this 
area when grown on 4 different major soil 
types. According to our plans we will have 
one nursery on each of the 4 different soil 
types which will be rather widely distributed 
in this area. We feel that it would facilitate 
our work somewhat if we would use the same 
random distribution in each of 4 locations. 
Why would it not be desirable to use the same 
random plan in each of the 4 locations? 


ANSWER: The probability values given in 
the standard tables for tests of significance 
were calculated from certain assumptions 
about the nature of experimental data. One 
of these assumptions is that the experimental 
errors which affect different units are inde- 
pendent of one another. This assumption may 
not be justified. In fact, in field experiments, 
the yields of neighboring plots are usually 
found to be positively correlated; which 
means that the experimental errors on neigh- 


boring plots also are correlated. Randomiza- 
tion effectively avoids bias from this source 
and permits the data to be treated as if errors 
were independent. 

Experiments have the common character- 
istic that, when they are repeated, the observed 
effects of the treatments vary from trial to 
trial. This variation introduces a degree of 
uncertainty into the interpretation of the re- 
sults. The results are said to be subject to 
experimental error. Whatever the source of 
the errors, replication of the experiment de- 
creases the error associated with the average 
effect of any treatment, provided that certain 
precautions are taken. These precautions 
must ensure that one treatment is no more 
likely to be favored in any replication than 
another, so that errors affecting any treatment 
tend to cancel out on the average as the num- 
ber of replications is increased. One essential 
safeguard is that the treatments be assigned to 
the experimental units at random. That is, 
the purpose of randomization is to avoid any 
systematic bias arising from the differences 
among the experimental units. Every treat- 
ment receives an equal chance of being as- 
signed to a favorable or unfavorable set of 
units. If treatment 1 and 4 appear adjacent in 
two replications in one location, and then the 
same randomization is used in four locations, 
you may have a biased pooled estimate of 
error for testing differences between these 
two treatments, although the estimate of the 
mean difference might be very accurate. 

The purpose of randomization is to guarantee 
the validity of the test of significance, this test 
being based on an estimate of error secured by 
replication. Therefore, it is not desirable to 
use the same random plan in each of 4 loca- 
tions. Gertrude M. Cox 


ANNUAL MEETING OF THE BIOMETRICS SECTION 


The annual business meeting of the Section 
was held at the Statler Hotel in Cleveland at 
12:30 on January 25, 1946, with 28 members 
and 13 guests present. Chairman Bliss report- 
ed on activities and development since the last 
previous annual meeting (September, 1944). 
This report appeared in the December issue of 
the Biometrics Bulletin. 


The nominating committee, consisting of 
Prof. W. G. Cochran and Miss Besse Day, re- 
ported their nominations. In the absence of 
additional nominations, the committee’s nom- 
inees were elected unanimously. 
of the Section for the new year are: Chair- 
man, D. B. DeLury; Secretary, H. W. Norton; 
Section Committee, E. J. deBeer, A. E. Brandt, 
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The officers — 
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J. W. Fertig, J. G. Osborne, and J. W. Tukey; 
Editorial Committee, Gertrude Cox, Editor; 
C. I. Bliss, W. G. Cochran, Churchill Eisenhart, 
F. R. Immer, H. W. Norton, G. W. Snedecor, 
and C, P. Winsor. 

Chairman Bliss outlined the plan for re- 
organization of the American Statistical As- 
sociation and opened discussion of similar 
plans for the Section. President Shewhart of 
the Association commented briefly, suggesting 
that the Section move in the direction of 
autonomy, the main advantage being indepen- 
dent control of its activities, especially in 
technical matters, but that it retain close ties 
with the American Statistical Association as 
the central organization of all statistical 
groups. Secretary Norton, speaking for A. E. 


Brandt, Chairman of the Constitution Com- 
mittee, presented a draft Constitution for the 
Biometrics Society, which it is proposed. that : 
the Section become, and indicated some of the 
questions which must be decided before a 
satisfactory constitution can be written. The 
discussion which followed centered around the | 
advantages and disadvantages of organizing as 
a Society rather than as a Section. It was 
emphasized that discussion of this question 
and of the constitution should be completed 
in 1946, so that the question can be voted upon 
at the next annual meeting, Following ap. 
proval of a resolution of commendation and 
appreciation of the work of the retiring chair- 
man, the meeting was adjourned. 


NEWS AND NOTES 


The Executive Committee of the Animal 
Vitamin Research Council (AVRC) met in New 
York City on December 15, 1945. The after- 
noon session was open to the general member- 
ship and included a report by C. I. BLISS, 
Chairman of the Statistical Committee on the 
analysis of data from a collaborative study by 
the AVRC of suggested rations of the A. O. 
A. C. chick assay for Vitamin D. ._ EL 
LeCLERG, pathologist, Bureau of Plant In- 
dustry, has been stationed at Baton Rouge, 
La. for several years. He is leaving the U. S. 
D. A. this spring to work in the Executive office 
of the President in the Bureau of Budget as a 
budget examiner. His duties will be those of 
_ a liaison officer between the research agencies 
of the U. S. D, A. and the Bureau of Budget. D. 
M. SEATH and Mr. LeClerg have been giving 
a course in statistics to a group of experiment 
station workers in the fields of nutrition and 
animal sciences. WILLIAM G. MAD- 
OW left the last of January for Sao Paulo, 
Brazil, to serve as Visiting Professor of Sta- 
tistics at the University of Sao Paulo for the 
full academic year which begins on March 15. 
He expects to return to the United States early 
in January, 1947, . S. C. SALMON and 
VICTOR R. BOSWELL have gone to Tokyo 
and other points in Japanese occupied terri- 
tory to help with problems relating to food 
production in the area. Mr. Salmon will deal 
with cereal and field crops while Mr. Boswell 
will be helping with fruit and vegetable crop 
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problems. A report will be expected when they 
return to the U. S. ROBERT PEN- 
QUITE is now with the Poultry Department 
at Iowa State College. . RALPH E. COM- 
STOCK, Plant Science Statistician, Institute of 
Statistics, recently spent ten days at Piedras, 
Puerto Rico, to review the designs and objec- 
tives of the Animal Husbandry research pro- 
jects of the Puerto Rico Agricultural Experi- 
ment Station. . CARL F. KOSSACK spent 
the first part of the war at the University of 
Oregon teaching under the ASTP program. 
During the last year he was in Washington, 
D. C., as a consultant with the Office of Field 
Service of the O.S.R.D. and was assigned to 
work for the Air Forces. He has returned to 
the University of Oregon. ERIC_ C. 
WOOD, Virol Ltd. (food specialists), Hanger 
Lane, Ealing, England, writes, “The Society of 
Public Analysts-and Other Analytical Chemists 
agreed in principle some time ago to the for- 
mation within itself of Groups for the study 
and furtherance of knowledge of such special- 
ized branches of analysis as might well be of 
interest only to a certain proportion of its 
members. A Biological Methods Group has 
recently been formed, having as its field the 
use of biological methods in chemical analy- 
ses.” This group held its inaugural meeting 
on October 17, 1945, at which time it was 
heartily agreed that their meetings would have 
a special interest to statisticians, since their 
work was so largely complementary to that of 


the bio-assayist. . . ALBERT LEE SCHRA- 
DER, Dept. of Horticulture, University of 
Maryland, College Park, Md., taught horti- 
culture at Shrivenham American University 
for awhile, then went to Amsterdam to discuss 
research activities in Agricultural Science, En- 
gineering and Commerce. He states, “I hope 
to get a number of new ideas of interest from 
the Dutch people.” During his stay in Eng- 
land, Mr. Schrader had the opportunity to 
visit several research stations. . . The work 
of D. BOYD SHANK, University of Arkansas, 
Fayetteville, consists of breeding and variety 
testing on corn and cotton. It is interesting 
to note that plant breeders now speak of the 
incomplete block designs as “conventional”. 

. WAYNE F. FREEMAN, Bureau of 
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Plant Industry, Experiment Station, Tifton, 
Ga., says the same thing: “My designs of ex- 
periments do not vary much from the now con- 
ventional lattices, triple lattices and lattice 
squares.” Mr. Freeman finished his Ph.D. at 
Illinois last spring. . . We were glad to 
hear from C. H. GOULDEN, Dominion Rust 
Research Laboratory, Winnipeg, Manitoba. 

. It has been an effort to see that some 
of the biological statisticians become sub- 
scribers to this Bulletin. That and secret war 
work has made it difficult to secure articles. 
Send us one, but remember, we want a discus- 
sion of the use of a statistical technique to 
biological research data. Inclusion of a nu- 
merical example is urged. 
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BULLETIN: Gertrude M. Cox, Chairman; C. 
I. Bliss, W. G. Cochran, F, R. Immer, J. Ney- 
man, H. W. Norton, L. J. Reed, G. W. Sned- 
ecor,, Sewall Wright. 

Material for the BULLETIN should be ad- — 
dressed to the Chairman of the Editorial Com- 
mittee, Institute of Statistics, North Carolina — 
State College, Raleigh, N. C., material for 
Queries should go to “uQeries,” Statistical 
Laboratory, Iowa State College, Ames, Iowa, — 
or to any member of the committee. 
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